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The values of the presently available truncated perturbative expressions for the pressure of the 
quark-gluon plasma at finite temperatures and finite chemical potential are trustworthy only at very 
£ — ■ large energies. When used down to temperatures close to the critical one (T c ), they suffer from large 

uncertainties due to the renormalization scale freedom. In order to reduce these uncertainties, we 
perform resummations of the pressure by applying two specific Pade-related approximants to the 
O^l ■ available perturbation series for the short-distance and for the long-distance contributions. In the 

' two contributions, we use two different renormalization scales which reflect different energy regions 

contributing to the different parts. Application of the obtained expressions at low temperatures 
^5 ■ is made possible by replacing the usual four-loop MS beta function for a s by its Borel-Pade re- 

summation, eliminating thus the unphysical Landau singularities of a s . The obtained results are 
' remarkably insensitive to the chosen renormalization scale and can be compared with lattice results 

t-H , - for the pressure p, the chemical potential contribution Ap to the pressure, and various susceptibil- 

ities. A good qualitative agreement with the lattice results is revealed down to temperatures close 
CO ' to To. 

> : 

O . PACS numbers: 12.38.Cy, ll.10.Wx, 12.38.Bx, 12.38.Mh 

m ■ 

<n : 

\Q • I. INTRODUCTION 

o ; 

The behavior of QCD matter at nonvanishing temperature and (quark or hadron) densities can be approached 
theoretically from two sides, both technically and kinematically. On the one hand, QCD-lattice calculations allow an 
analysis of the kinematic region of rather low temperatures T (around but above the phase-transition temperature 
T c , presently T < 5T C ) and even lower values of the chemical potentials (if (/ denotes the different quark flavors). 
■ The latter restriction is due to the fact that nonzero values of the chemical potential /i render the weight function in 
the partition function complex thus not permitting direct application of the standard Monte Carlo techniques. There 
are several tricks to circumvent this "sign-problem" : reweighting [l[ , Taylor expansion around /i = ||, analytic 
continuation from imaginary fi values 0], canonical formalism All of them are trustworthy only for small /i values. 
On the other hand, the behavior at high T and /i < lixT is expected to be described reliably by finite temperature 
(and density) perturbation theory (FTPT). Significant progress has been made within this latter approach during 
the last twenty years and the calculation of the thermodynamic potential (free energy or, equivalently, the pressure 
function p) has recently been pushed forward to the four-loop level, both for /i = [f| and for finite chemical potential 
Q. This is a big achievement because the corresponding truncated perturbation series (TPS) is in powers of the 
QCD-strong coupling parameter g rather than a = g 2 /(Air 2 ) = a s /ir, due to the well-known necessity of removing 
finite temperature infrared (IR) divergences by resumming the essential IR-sensitive diagrams ( "daisy diagrams" ) to 
all orders. The final series is up to 6th order in g, and that is essentially all that could be expected from perturbation 
theory since terms proportional to g 6 (not g e In g) and higher include genuine non-pcrturbative contributions which 
could be accessed only by non-perturbative methods, e.g. lattice calculations. 

Of course, any perturbatively obtained result is expected to represent the true physical situation only when the 
coupling parameter is small enough, which is the case in QCD at sufficiently high temperatures. Nevertheless, 
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evaluations of the 6th order result (for p and for Ap = — p(0)) down to temperatures below 1 GeV have been 
made in the literature, and the results could be brought in accordance with non-perturbative lattice data (available 
for T up to 1 GeV). At first sight this could be considered a triumph of finite temperature perturbation theory, and 
might tentatively be attributed to the higher order available. A second and more careful glimpse, however, reveals 
several sobering observations. The first problem is that the convergence behavior of the truncated perturbation series 
is manifestly weak: if increasingly higher orders in the series are added, the corresponding partial sums are changing 
wildly, jumping up and down. Only the step from 5th to 6th order has signs of moderation, but it is by no means clear 
whether this happens by chance or indicates a systematic improvement of the situation. Further, a common relatively 
high renormalization scale (RScl) 1 v w 2ttT (3> T) is used in such evaluations. For comparison with independent 
lattice data, the TPS results are used down to small T values where the lattice results are available. But here, another 
problem occurs: in each available order the corresponding TPS shows a dependence of the chosen RScl v, which is 
particularly strong for lower temperatures, thus making a consistent comparison of perturbative results (when used 
at such low energies) with lattice data doubtful. In fact we have no clear physical motivation for the best choice of v. 
In general, v should be chosen such that large (momentum dependent) logarithms in the TPS coefficients are avoided, 
which means - within an asymptotically free theory like QCD - that one is always on the safe side if v is taken to be 
near the lowest energy scales involved in the considered quantity. For the quark-gluon plasma at (high) temperature T 
(and for given chemical potentials /iy), what is the appropriate energy scale? Usually in the literature, the energy 2ttT 
of the lowest nonzero mode is taken as a measure since it determines the average energy of the constituents. But due 
to the collective effects additional (lower) physical scales are generated, namely the electric and magnetic screening 
masses (being of order gT and g 2 T, respectively). So what to choose for vl This question is not a purely theoretical 
one, but of considerable practical importance, because of the mentioned strong dependence of the TPS on v. And 
even if we neglect this intertwining of different scales within the thermodynamic quantities, and stick to 2irT as the 
relevant scale, we observe that a change of v by a single factor of 2 implies such a wide variation of the perturbative 
results at T ~ T c that no firm conclusion can be drawn as to the matching with lattice data (cf. Figs. [6] and [7]). It is 
clear, therefore, that no successful matching procedure can be obtained unless the strong RScl dependence is pinched 
down by improving the perturbative results. 

Within the present paper we offer a way for avoiding this unwanted (and unphysical) ambiguity by applying 
perturbation-theory-improving resummations of the basic TPS's. In this way we obtain what we consider a more 
reliable, but still perturbation-theory-based, description of the interesting thermodynamic quantities (here the pressure 
p) which, among other things, allow for a more credible comparison with lattice results. The method rests on replacing 
the (partially resummed) TPS by approximants which are much more stable under the variation of v than the TPS's 
themselves. It is well known that specifically Pade approximants Q for physical (measurable) quantities are stable 
under v variation Q, whereas similar improved Baker-Gammel approximants show this stability exactly [§]. 2 In a 
recent work [Tl| we have utilized Pade-related approximants to produce RScl-stable expressions for the pressure of the 
quark-gluon plasma at finite (large) temperature and zero chemical potential. Here, we apply the same technique to 
the case of finite nonzero quark chemical potential. In both cases, the gratifying fact is the high order of the available 
TPS which allows the use of higher-order Pade-related approximants and the choice of the most appropriate ones (see 
later). 

From the technical point of view we encounter a specific problem when applying Pade-related approximants directly 
to finite temperature perturbation theory. This is due to the fact that two different (infinite) classes of diagrams get 
involved in the whole mechanism: on the one hand those whose resummation is necessary for taming the finite- 
temperature IR divergences and which lead to contributions in powers of g (the so-called "daisy diagrams"); on 
the other hand the other diagrams, which give contributions in powers of g 2 and whose conversion from the TPS 
(polynomials) into Pade approximants (rational functions) had been known to result in RScl stability (which is 
exact in the large- /3q limit in the case of the diagonal Pade approximants [8]). Therefore, care has to be taken to 
avoid double counting when performing both resummations. As we have shown in Refs. a safe method consists 
in first decomposing the pressure into two parts, one containing the low-energy (effectively zero) modes and being 
responsible for the long-distance behavior of the correlation functions, and the other stemming from the nonzero modes 
and determining the short-range physics. Since both regions are, in principle, separately accessible by experiments, 
the corresponding expressions are physical in the sense that they should be independent of the renormalization scale 
v (at least up to the specified order in g). Further, no (infinite) resummation enters into the perturbation expressions 
of those parts when dimensional reduction is applied [l2| . Therefore, one can safely apply Pade-related resummation 



In our calculations we will denote the renormalization scale as v, keeping the symbol fj, for the chemical potential. The reader should 
be aware of the fact that in Refs. 0] and [f| the common RScl (in MS scheme) is denoted by JL and A, respectively. 
2 In addition, an extension of such RScl-independent approximants can be constructed, giving results which are simultaneously RScl- 
and scheme-independent [id . In the present work, the MS scheme is being used throughout. 
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to both parts independently and thus obtain two expressions which are both almost independent of the RScl choice. 
When using them down to low energies (T ~ T c ), we face another problem: the coupling parameter a(Q 2 ) acquires 
unphysical (Landau)singularities at low energies Q < 1 GeV if using MS TPS f3 function. We circumvent this problem 
by using again appropriately resummed versions of (3. We finally end up with expressions for the pressure and for other 
measurable quantities (quark number susceptibilities) which are almost free of the RScl uncertainty and therefore apt 
for comparison with lattice data. 

In Sec. II we present the perturbative results on which our analysis rests and describe how to perform the physical 
separation of the pressure into the long-range and the short-range part. Section III contains an analysis of the 
possible resummation procedures which leads to the optimal choice. We then present the numerical results, where 
we put the main emphasis both on the effect of finite fi values and on the comparison with the corresponding lattice 
data. Appendix lAl compiles basic formulae and expressions, available in the literature and adapted to our approach, 
including expressions for the coefficients of the TPS's. Appendix [Bl describes the method which allows us to extrapolate 
the QCD renormalization group equation to sufficiently low energies, thereby circumventing the unphysical Landau 
singularities of the coupling parameter. 



II. SEPARATION OF LONG- AND SHORT-DISTANCE PRESSURE 



Our starting point are the FTPT results for the pressure of the quark-gluon plasma which, for a homogeneous 
system, is equal to the free energy per volume (up to a sign difference) - considered as a function of its temperature 
T and the chemicalpotentials of the various quark flavors fif (/ = 1, ...n). It has been calculated up to 0(g 6 lng) by 
K. Kajantie el al. [5[ for the case of vanishing chemical potentials, and for the general case (/if ^ 0) by A. Vuorinen 
Q . These high order results include the summation of an infinite class of certain diagrams - necessary for taming an 
IR singularity which occurs only at T ^ 0. The final results could be practically achieved only by a technical trick, 
namely by separating the energy-momentum region of the contributing modes into three parts, characterized by the 
momentum scales 27rT, gT and g 2 T, such that the full pressure p is decomposed according to 

P = Pe + Pm + Pg ■ (1) 

Note that this decomposition makes strict sense only for high enough temperatures where g(T) < 1. Here, pe 
represents the contributions of all degrees of freedom associated with the nonvanishing Matsubara modes, whereas 
Pm + Pg comprises the contributions of the zero modes (of bosonic fields), thereby implicitly representing also the 
necessary sum over all (daisy) diagrams. The latter ones are static modes, hence their contributions can be effectively 
described by a three-dimensional (in general d-dimensional) purely bosonic field theory (dimensional reduction [12[) 
determined by the electrostatic QCD (EQCD) Lagrangian 

£eqcd = \^F 2 + Tr[A, io] 2 + m 2 E TvA 2 + X^(TrA 2 ) 2 + X^TrA* + Z J£_ ^ ^ Trig + ... (2) 

Here Aq denotes an effective (<i-dimensional) scalar field and the Ai (i = 1, ...d) define a d-dimensional vector field, 
both in matrix notation (i M = A^T a ); 3 Di = di — ig-^Ai] Fij = (i/gE)[Di, Dj]. The parameters of this effective theory 
are the (electrostatic) screening mass toe (~ gT), the effective coupling parameter g 2 , (~ g 2 T) and the four- vertex 

couplings g A T). In the case of d = 3, A E and are not independent and one can choose = 

- this will be done in the following. There are additional coupling parameters connected with Lagrangian operators 
of higher dimensions, since £eqcd defines a non-renormalizable theory which makes sense only for momenta below a 
certain (UV)cutoff Ae- In our case Ae separates the region of momenta ~ 2-7rT from the momenta ~ gT and smaller. 
The effective parameters can be connected to the parameters of the underlying QCD by means of the well-known 



3 The fields Ag, A? entering here are not identical to the gluon fields in £qcd but are the effective fields obtained after the high-energy 
modes have been integrated out. 
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matching procedure 13] yielding 



= T 2 l g 2 



A, 



A^ In ■ 



2vrT 



-4. 



)+0(e 2 ) 



(4tt) 2 



Tig 2 



(4tt) 2 



[/3 E 4 + 0(e)]+0(<? 6 



A® = 



(3) 
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Here, e = (3 — d)/2; coefficients A4-A7, A^-A^ and /3e4 are complicated functions of the chemical potentials, the 
latter appearing in the coefficients via the dimensionless quantities Jit = fif /(2ttT). Their expressions are collected 
in Appendix [XI together with other coefficients to appear in Eqs. (f5|)- (fTU|) . Note that the common RScl v c appears in 
these expressions. The effective mass Toe ~ gT arises due to the color-electric screening. 

Since there is, in addition, color-magnetic screening at energies proportional to the corresponding magnetic screening 
mass Tom ~ 9 2 T, the long-distance part of the pressure can be further subdivided into pm and pc, where pm is 
determined by £eqcd and pa by the (magnetostatic) Lagrangian 



£mqcd — 2^ r ^' 



(6) 



with Fij — [i/gu)[Di, Dj] and Di — di — igyiAi. This Lagrangian defines the effective theory for energies below Am, 
i.e., for energies ~ g 2 T and smaller. A similar matching procedure as before determines <?m in terms of the parameters 
of the higher-energy Lagrangian £eqcd and gives 



g 2 M = g 2 K + o(g 3 ). 



(7) 



In the case of nonzero chemical potentials, two scales get involved (T and p, = p,f) and, therefore, the concept of 
dimensional reduction is expected to be applicable only if the magnitude of the chemical potentials is small compared 
to 2nT [lij . From comparison with numerical results for correlation lengths, it is expected that the restriction /1 < 4T 
is safe. 

Based on Lagrangians @ and ©, and on the ordinary QCD-Lagrangian, the various parts of the pressure have been 
calculated perturbatively by Vuorinen [fa. The calculations are based on dimensional regularization with a common 
RScl v c (the notation A is used in Ref. Q for the common RScl). 

The result for p& is 



PE 

T 



A!+g 2 (A 2 + 0(e)) 



(4tt) 2 



1 



U \ -6A 4 + 4^1n^ + A 3 + 0(6 



(47T) 



: g 6 (fci+0(e))+0(g 8 ) 



(8) 



Coefficients Aj (i = 1,2,3) and A^ are collected in Appendix [XJ Coefficient /3ei at g 6 is still unknown. However, 
/?ei must include a term proportional to 1/e. The 1/e-terms will be disposed of in the following because such terms 
must cancel in the sum (JT|), and the finite part of /3ei will contain a free (adjustable) parameter later in this work. 
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The results for pm and pg, which can be obtained from the effective Lagrangians Q and ([6]), respectively, are 
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Here, tom is the magnetic screening mass tom = ^9m (~ 9 2 T), and the notation ~ftf = fj,f/(2irT) is used. We 
observe that the chemical potentials show up both explicitly via JLj and implicitly via g E and m E . Coefficients 
ttMi, ckM2, olq, (h/ii, 0M2 are independent of fif (see Appendix |A"|) . and only /3m2 is not yet known. In Eqs. ([9]) and ([101) 
we have already separated all "divergent" terms (proportional to 1/e) from the finite contributions. The quantity 
Pm/T starts effectively with order g 3 and pg/T with order g e . The full ~ g e term in pq/T cannot be determined in 
a perturbative way; however, Sq was estimated in Ref. [l5[ to be i5g = —0.2 ± 0.6. Here, it will be treated as a free 
parameter within the limits Sq — ± 1. On the left-hand sides of Eqs. (|8))- (fT0)) . the common denominator T must in 
fact be replaced by Tz^~ 2c ; however, the common factor v~ 2e reduces to unity in the e — * limit, and will be ignored 
because the sum p E + Pg + Pm is finite in this limit. 

Our main task now consists in deducing from these formulae expressions for the physical short-range and long-range 
parts of the pressure. Thereby the word "physical" indicates that they lead to measurable effects (for instance, the 
long and/or short-range behavior of static correlation functions, etc.). We proceed in several steps (see also Ref. [HI]): 

(i) Regularization: 

The long-range part of the pressure is represented by pm+Pg = Pm+g (both are due to zero mode contributions). 
We regularize pm+g by adding to pm+g/T the following counterterm: 



CT= -Jo 
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(11) 



The same counterterm has to be subtracted from Pe/T. By expanding the effective-theory-parameters g E and 
m E in powers of g [cf. Eqs. Q, Q], one can explicitly show that in this way the 1/e-term of order g 4 in 
expression ([5]) for p E /T gets cancelled. Further, the 1/e-term of order g 6 included in the otherwise unknown 
Pei coefficient must also get cancelled when counterterm (fTTj) is subtracted from pe/T of Eq. ([5]). Counterterm 
(|lip contains one finite term of 0(g A ) stemming from (l/e)g E m E since m E includes a term proportional to g 2 e, 
Eq. ([3]). This finite term then shows up in the new subtracted expression for pe/T. In addition, counterterm 
(llip contains several finite terms of 0(g e ) which show up in the new subtracted expression for p E /T. Finally, 
the limit e — > can be performed, yielding finite results both for p E and for pm+g- 

(ii) Reconstruction of the factorization scale and introduction of various RScl's: 

We have already noted that a common renormalization scale (denoted here as v c ) has been used for the per- 
turbative calculations leading to Eqs. (f5|- (110p . On the other hand, when constructing the physical long- and 
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short-range contributions, respectively, we decompose the whole energy range in the way addressed before, 



Therefore, these factorization scales Ae and Am which separate the different energy regions have to emerge in 
the physical expressions - but only in such a way that, when adding all three contributions they completely 
disappear. We bring the Ae factorization scale to light in the following way: the scale v c in expressions ([8])- 
(|10jl . where pe/T and pm+g/T are modified in the afore-mentioned way by the counterterm (jllj) . is interpreted 
simultaneously as the factorization scale Ae and as the common RScl v c . We then evolve g(y = Ae) in pe 
to (/(^e) where v — ve is a new, physically more adequate higher RScl: z^e ~ 27rT. On the other hand, in 
Pm+G: we evolve g(y = Ae), which appears implicitly there (explicitly in the RScl-independent gE, mE, Ae"*), 
to g{vM) where v — z/m is a new, physically more adequate lower RScl: vm ~ We (~ gT). The evolution is 
performed according to perturbative renormalization group equation (RGE), requiring RScl independence of 
Pe, on the one hand, and of cz E and m E (and thus of pm+g), on the other hand, since all these quantities are 
physical. In pe this results in In z^E-dependent terms in the coefficients of TPS, and in g 2 , and m E (which enter 
Pm+g) this results in In z^M-dependent terms in the coefficients of their TPS's. The coefficients of expansion 
of pe in powers of g 2 (vE) then have explicitly the "genuinely" AE-dependent parts, and RScl-dependent parts. 
The coefficients of expansion of pm+g i n powers of (gg/mg) remain unchanged, with Ae dependence as before, 
while the coefficients of expansion of m E and g| in powers of g 2 (vM) obtain RScl dependence and, at the order 
considered, lose Ae dependence. 4 Formally, the two RScl's ve and um can take on arbitrary values in these 
expressions. 

(iii) Determination of the AE-dependent part of coefficient at 0(g 6 ) in pe'- 

The sum pe + Pm+g has to be independent of the factorization scale Ae- It can be checked explicitly that 
this independence is true up to 0(g 4 ). Although the coefficient at g 6 in pe is not known, its Ae dependence is 
dictated by the condition of Ae independence of pe + Pm+g at 0(g e ), when the latter quantity is expanded in 
powers of a common g = g(v). Note that Ae dependence of pm+g at 0(g 6 ) is known, cf. Eqs. (|5|)- (fTUl) . Further, 
v dependence of the coefficient at g 6 in pe is known from the requirement of the RScl independence of pe- The 
remaining unknown part of the coefficient at 0(g e ) in pe (we will call it Se) is then independent of Ae and of 
v, and can again be freely adjusted. 

Performing these steps we finally obtain the following form of the physical decomposition of the pressure function into 
short- and long-distance parts (we denote these physical parts by a bar) 



The long-distance part, representing the contributions of momenta below the factorization scale Ae (when gT < Ae < 
2ttT) takes the form 



namely 



0(g 2 T) < A M < O(gT) < A E < <D(2nT) . 



(12) 



P = PE +PM+G- 



(13) 




(14) 



4 This is not so in the scalar g 2 </!> 4 theory, where the requirement of RScl independence of the Debye screening mass yields, at 0(g 4 ), a 
residual dependence on the factorization scale (cf. Ref. second entry). 
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where we used notation (|A1|) - (| A2|) for the chemical potential parameter The parameters toe,3e and A E ^ of the 
effective theory EQCD are defined by their expansion into powers of g[v) = g 



fAj 1 + (J-) 2 Ata ( 
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Here, /3o = (1/4)(11 — 2n//3) is the one-loop QCD RGE coefficient, n/ being the number of active quark flavors. 
The RScl v in Eqs. (fl"5|) - (fTT|) will be v = vyi (~ hie ~ gT), and we fix it according to relation Q27p. Note that, 
by the afore-described procedure, the expansion coefficients of to e and <7 E at the considered order do not have any 
dependence of the factorization scale Ae, but are RScl-dependent. The last term in expansion (fT4]) involves the 
third EQCD matching parameter A^ p^. Il7|. which is independent of the chemical potentials This term can be 

expressed as a power series in powers of g, but only the leading term is known A^/toe oc g 3 . We prefer to express it 
in powers of the EQCD parameter <7 E /m E (Ref. [ll|) 



(l)_2 1 
"3(4tt)2 



(9-n f )A 4 m E ( 

1 TOE 




(18) 



While coefficients aMi, etc, Pmi in Eq. (fT4)) and Ai (i = 4, ...7) in Eqs. (fl"5)) - (fl"6|) are known and collected in Appendix^] 
parameter Sq is well estimated [15| . and /?M2 is unknown. The dependence on the chemical potentials fj,r (/ = 1, n) 
in Pm+g appears explicitly in the term proportional to Jl\ and implicitly via the parameters toe,<?e which contain 
/^/-dependent coefficients A$. The RScl v — z^m appears in Pm+g implicitly, via expansions p5[) and (| 1 6[) for m| and 

ah 

The physical short-distance part, determined by the energy- momentum range above the factorization scale Ae, can 
be written in a dimensionless form as follows: 



1 

y>4 



p E =A 1 +4n 2 A 2 R E an , 



(19) 



where i? E an denotes the canonically normalized perturbation series in powers of g = g{v): 
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Here, ft = (1/16) (102 - 38n//3) is the two-loop RGE coefficent, and v 
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k = 27rexp 
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2wT) is the short-range RScl, and 

(21) 



has been introduced for obtaining compact expressions. The ^/-dependent function IC3 has been obtained by requiring 
that the dependence of p~m+g on the factorization scale Ae at order g e cancel with that in p E . It arises in the following 
way: inserting expansions (fl"B")) - (fTT)) into (fT4")) yields to order g 6 the As-dependent term 
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^4 Pm+g (5 hiA E ) = 
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Cancellation with the corresponding term in pe is obtained therefore if 



-(A 6 + A 4 A 7 ) - 18(omi + a G ) + 7^ n Wi 
! 27 



(22) 



(23) 
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The remaining (unknown) coefficient at g e within pe, which we denote by Se, is independent of Ae; it may, however, 
depend on the (small) ratios ~pt — fif/(2nT). 

Since the true values of p~E, 3 E , m E and pm+g are RScl-independent, this motivates us to apply Pade-related 
resummations separately to the TPS's of these quantities, yielding expressions which are much less RScl-dependcnt 
than the corresponding TPS's. In this context, we note that the RScl independence of perturbation expansions (with 
infinitely many terms) for the afore-mentioned quantities is only formal, because these series diverge. They diverge 
strongly at low temperatures when g{v^) [and even more so g{vyi)\ gets large. Thus, at low temperatures, the formal 
RScl dependence of the series does not help in direct evaluations of the TPS's since the (higher-order-)RScl-dependent 
corrections to the TPS's are large. Therefore, for this kinematical region, the only way out of this RScl dependence 
dilemma seems to be the conversion of the perturbative expressions to other approximants which are much less RScl- 
dependent at each finite order, even at low temperatures. And exactly that is the main motivation for applying 
Pade(-related) resummations. 

Before doing so, we discuss the other uncertainties of the TPS's, namely the uncalculated constants /3m2 and 5e, 
and try to estimate their expected size. The constant /3m2 appears in expression (I14p in combination with the term 
I^Ae/toe), where the latter is expected to dominate. Therefore, |/3m2| < 15 re pre sents a rather generous uncertainty 
bound for (3m2- On the other hand, the constant Sq was estimated in Ref. [15| to be <5g = —0.2 ± 0.6. Here, it 
will be treated as a free parameter within the limits Sq = Oil. Thus, we will allow the following variation of the 
aforementioned parameters: 

- 1< 8 G < +1 , -15 < (3 M 2 < +15 . (24) 

Concerning Se [see Eq. (|20[) ]. we note that the parameter k was introduced in such a way that the hi[is / (2ttT)]- 
independent part of the coefficient at g in p~E is absorbed by a term proportional to ln[AE/(KT)]. The coefficient at 
g was then organized into a polynomial in powers of the aforementioned two logarithms. It is reasonable to expect 
that the ln[^/(27rT)]-independent part of this coefficient is absorbed to a large degree by a term proportional to 
ln[A E /(KT)]. Therefore, parameter <5e is expected to be small and the following variation of this unknown parameter 
appears to be rather generous: 

- M < S E < +N , (25) 

where 

^ 6 ^ 3 !„(§). (26) 

Parameter 5e depends only on the small parameters Jlj. On the other hand, the bounds ±|fc2| will have an additional 
slight dependence on temperature T because we take Ae = \j2-KTmE(T) (~ g 1 / 2 T). 

Formulas in this Section are in close analogy with those in our previous work [llT | , invol ving now additional (small) 
parameters Jlj = n f / (2ttT) . Furthermore, up to terms g 4 , they coincide with those of Ref. |X3j] when Jjj = 0. Further, 
rccxpanding Pm+g °f Eq- 03) m powers of the coupling parameter g = g{v), and adding it to expansion (TT5)) - (f2T))) for 
p E while using there the same RScl v, gives the same expansion as the one obtained in Refs. Q for pe+m+g- 

III. RESUMMATION AND NUMERICAL RESULTS 

Our next step is to apply specific Pade related resummations to evaluate separately the long-distance (|T4])- (jT8j) and 
short-distance (|H)]) - (|20j) contributions to the pressure. 

In principle, we could utilize Pade (P[N/M]) or Pade-Borel (PB [N/M]) approximants 5 of any possible order [N/M] 
which is compatible with the order n (the highest power of expansion parameter) of the TPS: N + M < n. So we have 
a certain freedom of choice. We use it for achieving physically desirable features. These are: 

(a) Significantly suppressed RScl dependence of both resummed pm+g and Pe, with the two RScl's z^m and ve 
varying in the regimes vm ~ gT and ve ~ 2?tT. Minimal RScl dependence is achieved in general for diagonal or 
near-diagonal approximants (N « M), thus we expect that such approximants will be preferred. 



5 For a short description of Pade and Borel-Pade approximants, see Appendix of Ref. fTljl . 
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(b) p E + Pm+g should have as little dependence on the factorization scale Ae as possible. Note that the sum of the 
original TPS's Eqs. (fH?)) - P0|) and ([H]) , when expanded in powers of a common g — g{v) up to ~ g 6 , is completely 
stable under variation of Ae- On the other hand, the individual parts show significant Ae dependence. Since 
these individual dependences get (individually) changed by resummations, we have to optimize the approximants 
in the sense of maximally reducing the artificial Ae dependence of p~E + Pm+g- 

(c) Pe + Pu+g should not surpass the value of Pidcai (pressure of the ideal gas), even at low temperatures close to 
the critical temperature T c (see Ref. [TTJ] for arguments in this direction). 

Since Pm+g is expanded in powers of (<? E /m E ), we first have to calculate the EQCD parameters m E and <? E . We 
evaluate them as Pade resummations P[l/l](g 2 ) of TPS's (fTSf and (fl"6|) . thereby banking upon the better convergence 
behavior of Pade approximants and reducing the unphysical RScl dependence dramatically. 

Within our previous paper , we have shown for the case of zero chemical potentials that, among the resummations 
of the Pade (P) and Borel-Pade-type (BP) of the perturbation expansions (fl4|) for pm+g arid (|20|) for i? E an , the only 
physically acceptable ones in the afore-mentioned sense (a)-(c) are: 

(I) the Pade approximant P[0/3] in terms of expansion variable g B /m,E for pm+g/^| 01 ^q- without the 
Ag-'-term. Note that, since g E / TO E ~ 9 and m E ~ g 3 , this emulates partly the diagonal P[3/3](gr) for Pm+g- Further, 
the Ag^-term in Pm+g/ to e 01 ^ c l- $M ^ s evaluated according to Eq. (JTSJ) , making it less RScl-dependent. 

(II) the Borel-Pade approximant BP[l/2] as a function of a = g 2 /(2n) 2 for R E an of Eq. $20)) (=> p E ). 

Explicit construction of the Pade's P[l/l](a) for m| and g E , P[0/3](g E /wE) for pm+g/^e' anc ^ BP [1/2] (o) for 
i? E an can be read off from Appendix in Ref. [HI]. 6 The A E ^-term is not included in the aforementioned Pade-related 
resummations, but is evaluated separately [Eq. (fl"8"]) ] and added, because it is expected to represent diagrams with 
new, different, topologies. 

Since the chemical potentials in expressions (114|) and (fT9l) -(|20 |) are assumed to be limited in the sense that ~pf < 1, 
they can be regarded as perturbations to the /I/ = case [6J. Therefore, we will apply the same approximants as in 
the p/ = case of Ref. i- e -> those mentioned above. Furthermore, just as in Ref. [ll|, we will fix the two RScl's 
pe and vm according to relations 



ve 



2ttT , ul A = ml{T,-p ]VM ) \= m^ ' 2 (T, /!)] , (27) 



where to e is taken to be the Pade approximant P[l/l](a) of expansion (|15|) . as mentioned before, with g = g(vy\). 
We note that now the long-distance RScl vm will depend on both the temperature T and the chemical potential 
H = 2irT~p. At sufficiently high temperatures, we have i>m ~ gT, cf. Eq. (|15p . The factorization scale A E is chosen to 
be just in-between the two RScl's (|27|) on the log scale 



A E = = V 2nTm { ° ) . (28) 

With all these quantities fixed, there still is a problem of obtaining results at low temperatures close to the critical 
temperature T c ps 0.2 GeV, where the values of the long-distance RScl vm are much below 1 GeV. At such scales, the 
usual perturbative MS couplant a{v 2 ) = [g(v)/(2ir)] 2 diverges as a result of the unphysical Landau singularities, the 
latter being the consequence of the beta function (3(a) occuring in the form of the (four-loop) truncated perturbation 
series (TPS). A partial remedy to the related problem of unreliability of evolution of a(y 2 ) at low v was presented in 
our previous works pT| , where we used Pade P[2/3](a) for (3(a). However, the problem of the Landau singularities 
in the low-energy space-like regime persists. Therefore, in Appendix [B] we present another resummation of the MS 
four-loop beta function, of the Borel-Pade-type (BP). We show that BP[2/2] and BP[l/3], in MS, result in evolution 
which keeps a(is 2 ) finite down to v 2 = 0. The two BP's give mutually similar results. Even more so, when varying 
the scheme, e.g., by changing the values of (32 and (3% coefficients by about 50%, the main qualitative features of the 
low-RScl evolution survive. We will adopt for a s (v 2 ,MS) = ira(v 2 ) (with the space-like RScl values q 2 = —v 2 ) the 
reference value 

g 2 (v — m T )/(47r) = u s (v 2 = m 2 ,MS) = 0.334 (n f = 2 or 3) , (29) 



Fig. 2(b) of Ref. [ill , which shows RScl dependence of various Pade-related resummations of p^,, has numerical errors for RScl values 
i^E 7^ 2ttT due to a mistake in one of our programs; corrected curves show that, in addition, P[l/2] resummation for R^ n has an 
acceptably suppressed RScl dependence. However, then p/pidcal > 1 at T ~ 1 GeV, making it unacceptable (where is from P[l/2] 
of R°*°-, and p M +G from P[0/3] of p M+G /m|). The conclusion in Ref. 11] that only BP[l/2] for #£f n and P[0/3] for p M +G/"i| are 
acceptable remains unaffected. 
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which is approximately the value extracted from the hadronic r decay data 0, [T^|; and we will evolve a(v 2 ) by the 
BP [2/2] beta function (see Appendix [Bl Figs.[T3Ta) and fTTT a) for n* = 3, 2, respectively). Sometimes, for comparison, 
the evolution by the BP[l/3] beta function will be used (Figs. [T37 b) and fTTT b) for rif = 3, 2). Details of the definition 
of BP[i/j] for the beta function and other details are given in Appendix IBl 

Most of the following calculations are performed for the case of two active massless quark flavors (rif = 2), in order 
to facilitate comparison with the lattice calculations of Refs. [2{J|2l|. Some of the calculations will be performed for 
rif = 3 in order to see the rif dependence of the results. We adopt the notations used in Ref. [21!]: 

Ap = p(T;fXf)-p{T;fx f = 0) , (30) 
N = oG"« + Md) , VI = o(Mu - Mrf) . (31) 



2 Kt ^ ' 1 ^ 2 
dp 



(32) 



d 2 p _ d 2 p 

Xq ~ d(fi q ) 2 ' XI ~ d{^y ' 



where the partial derivatives with respect to fix {X — q, I) are taken at constant T and constant \j,y (Y — I, q, resp.). 
Here, n q is the quark number density (at fij — 0); \q an d Xi are the quark number and isovector susceptibilities. 
When rif — 3, we will take fi s = 0, as in the lattice calculations of Refs. [22l. [23|. Further, we will use for the critical 
temperature the value T c = 0.17 GeV as in Refs. [20l [2l| . both for rif = 2 and rif — 3. Unless otherwise stated, 
the unknown parameters Sq, /3m2 [Eq. (|24|) ] and 5e [Eq. ([25]) ] will be set equal to their central value zero, and the 
RScl's fE and vm for the short- and long-distance parts of the pressure will take on the 'canonical' values according 
to Eq. (|2"T|) . Numerical calculations were performed using Mathematica [24[ ■ 

In our calculation of Ap we will resum separately, in the aforementioned way, p(T,fif) and p(T, fif = 0), and then 
subtract the two quantities. We prefer this approach (instead of trying various resummations of the perturbation 
series of the quantity Ap) because the directly measured physical quantities are the full pressures. 
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FIG. 1: Pressure p (divided by T 4 ) as a function of temperature, at various values for the ratios piq/T involving the chemical potential 
Hq when (a) Borel-Pade BP[2/2] and (b) BP [1/3] is used for the beta function /3(g%). 



Numerical resummations, performed in the way described above, give us for the total pressure p = Pe + Pm+g the 
results presented in Figs. [U for various values of the ratio jiq/T (where \iq — fi u = fi d ), as a function of temperature 
in the vicinity of T c . Comparison of Figs. [T^a) and (b) further reveals that the results do not change significantly 
when the type of the BP resummation of the beta function is changed. For better visualization, we present in Fig. [2] 
a three-dimensional image, showing p/T as a function of T/T c and of ii q fT c (fij = 0), for the choice of parameters, 
RScl's, and resummation approximants equal to that of Fig. [lja). Note, however, that in Fig. [2] the second axis is 
Hq/T c and not fi q /T (the latter quantity is kept fixed in the separate curves of Fig. [TJa)). 

In Fig. [3] we present the corresponding results for the pressure difference Ap = p(T;ii q ) — p(T;ii q — 0), for five 
different values of \i q jT (= 0.2,0.4,0.6,0.8,1.). We included, for comparison, the results of the evaluation of the 
simple truncated perturbation series (TPS) in powers of g[y — 2ttT) as dotted lines there. These were obtained by 
using for pE the TPS in powers of g 2 (v = 2ttT), and for Pm+g the TPS in powers of g(v = 2wT). The latter TPS is 
obtained by using expansions (|15p - (ll7p in powers of g(v) in expansion l|14p for Pm+g (also in the logarithms there), 
and setting the RScl v — 2ttT. The unknown parameters 5e, Sq, and /3m2, which affect these TPS's at 0(g 6 ), were 
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FIG. 3: As in Figs.[T] but now for the pressure difference Ap divided by T (solid lines). In addition, the TPS results are included as 
dotted lines; and the lattice calculation results of Ref. [2l| are included as crosses, where the depicted error bars include only specific 
statistical errors. 



all set equal to zero here. In addition, the TPS's truncated at 0(g 5 ) are presented, for the aforementioned five values 
of Hq/T. We note that such types of TPS evaluation (with the common high RScl ~ 2ttT) have been often used in 
the literature to evaluate p and/or Ap. The TPS's presented here do not diverge when approaching even very low 
values of temperature because for the beta functions we use BP [2/2] (or: BP [1/3] in Fig. [1Tb)). Fig. Oshows that our 
Pade-related evaluations, while being somewhat higher, reproduce the lattice results for Ap/T 4 to within 20%, even 
at low temperatures T « T c . On the other hand, the TPS results are very unstable under the change of the truncation 
order. It appears to be a coincidence that the TPS's truncated at 0(g 5 ) are in good agreement with the lattice data. 
Incidentally, the latter TPS's have values similar to those of C(g 6 )-TPS's with Se = [in the latter case, Pe has 
the coefficient at g 6 equal zero, cf. Eq. (p?5|) and (p?0"]) ]. Stated otherwise, the 0(g 6 ) TPS's are quite unstable under 
the variation of the unknown parameter Se while our Pade-related resummations are quite stable (see also Figs. [4] and 
[5]). Furthermore, the TPS results show a strong RScl dependence, as will be shown shortly. 

A general remark on the lattice data (included in Fig. [3]) and their significance is in order here. The quoted error 
bars denote only specific statistical errors and do not represent further uncertainties. The aforementioned lattice data 
have various uncertainties, among them the generic uncertainties of lattice calculations coming from the continuum 
limit effects (of up to 10%, cf. 20]), from the finite size effects (of about 5%, cf. [25j|). and from the uncertainties 
of the value of T c (of 2-3%). The results for Ap for finite chemical potentials suffer from additional problems: since 
finite /i g -values are treated in Ref. [2l[ by applying a Taylor expansion in powers of n q /T (cf. Eq. (3.1) of Ref. [2l|). 
one has to calculate the corresponding coefficients c„(T) - more of them when fiq/T is higher. In Ref. [2l[ only the 
first three non-vanishing coefficients (n — 2,4,6) have been calculated - with high instabilities already for n = 6 
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(cf. Fig. 1 in [21]). The error bars in the lattice data in Fig. [3l as well as in Figs. [8]|9l present only these specific 
uncertainties in calculation of the three c n 's. However, the terms with n — 8, 10, . . . are not included. Consequently, 
the lattice results for large fJ. q /T values (fJ. q /T > 0.8) have to be considered with reservation. An educated guess 
leads to the expectation that the data in Fig. [3] have an overall uncertainty of around 15%, and probably even higher 
when [iqjT w 1. Therefore, it is fair to say that our predictions are in reasonable agreement with lattice data down 
to T = T c , at least as long as [i q < T. Analogous statements are valid for the comparison with lattice data in 
Figs. [5] and 03 (see later). What we do not yet understand is the apparent systematics of the lattice results: they lie 
systematically below our predictions, in particular for high /i q values. Whether this demonstrates a lattice artefact, 
possibly connected with the rather large bare quark mass used there, has to be further investigated. In this context, 
we further note that the difference between our and lattice results could not be significantly reduced by choosing 
different values for the unknown parameters £e, /3m2 5 <>g (which were set equal to zero in Fig. [3]), at least if varying 
them within the generous ranges specified in Eqs. (f2"4"|) and (f2"5)) . In fact, variation of these parameters can decrease 
our results at 1 < T/T c < 2 by less than 1% (see Fig. [5] later). 
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FIG. 4: (a) Pressure p (divided by T 4 ) as a function of temperature, when the unknown parameters <5g, 0M2 and 5e are varied according 
to Eqs. i ]24[ ) and (1251 . for two different values of the ratio fiq/T: 0.8 and 0; (b) same as in (a), but for TPS's (with the common RScl 
v = 2-kT). 



In Fig. (Ha) we present variations of our results for p/T 4 in a wide temperature regime when the unknown parameters 
Pm2 and <5e are varied according to Eqs. (f2"4")) and ([23| . for two different fixed values of the ratio fi q /T (= 0.8,0.). 
In Fig. 0Jb), the analogous results for the aforementioned simple TPS's are shown. We see that our results are 
remarkably stable under the rather generous variations of the three unknown parameters, whereas this is definitely 
not the case with the TPS's. The dependence on the unknown parameters <5e and <5q is strong in the TPS's, while 
the Pade-related resummation results are almost independent of them. The dependence on the parameter /3m2 is too 
weak to be seen. 
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FIG. 5: As in Figs.gl but now for Ap/T 4 . 
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In Figs. 02a), (b) we present Ap/T 4 in the way completely analogous to the presentation of p/T 4 in Figs. IHa),(b). 
The conclusions for the (in)stability of the calculated Ap under the variation of the unknown parameters are similar 
to those for p. The independence of the parameter Sq in the TPS's in Fig. |4jb) is a direct consequence of the \i q 
independence of Sq. The dependence on the parameter /3m2 is weak in the TPS's, and too small to be seen in the 
Pade-type resummation. However, while there is almost no dependence on the unknown parameter <5e in the Pade- 
related resummation results, the dependence on <5e is quite drastic in the TPS's. This behavior also makes plausible 
the fact that, by adjusting the value of the unknown parameter Se, we can, in a way, fine-tune the TPS results to 
come close to the lattice results (see also Fig. |3J). 
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FIG. 6: Pressure p/T 4 as a function of temperature, when the renormalization scales and vyi are varied by factor two around the 
canonical values H27I I: (a) when rij = 2; (b) when nj = 3. The values of the ratios p q /T are either 0.8 or zero. Included are also the 
corresponding sets of curves with the TPS evaluation. 



After having shown that our results are fairly insensitive to the still existing unknown parts of the perturbation 
series (at ~ g 6 ), we now come to the most important results of our approach: the stability under variation of the 
(two) RScl's, even at very low temperatures. This is manifested in Figs. [5] and [71 

In Figs.[6]Ja),(b), we present the behavior of our evaluated results for p/T 4 when the RScl's z^e and Vyi are varied by 
factor two around the canonical values P7|) . for nj — 2,3, respectively. Two sets of curves are given, for \i q /T = 0.8 
and zero, respectively. In addition, the corresponding sets of curves for the aforementioned simple TPS's are shown, 
where now the common RScl v is varied from ttT to AnT. We see that our evaluated results for p are much more 
stable under the variation of RScl than the TPS results, down to very low temperatures T « T c . This is the same 
conclusion as the one obtained in our previous work for the case of zero chemical potential (fi q = 0) . For additional 
comparisons, we included in Figs.[5]the results of our Padc-rclatcd evaluation of p/T 4 with canonical RScl's (|2"T|) when 
the MS beta function /3(a) is (four-loop) TPS, and when it is Pade P[2/3](a). We see that the results in such cases, 
when they exist, almost coincide with the solid-line curves, i.e., with those with (3{a) = BP [2/2] (a). However, due to 
the Landau singularities of a(v 2 ) at low RScl's in the aforementioned cases of (3 = TPS or P[2/3] (cf. Appendix |B|) . 
the corresponding curves exist (i.e., do not blow up) only down to T m i n w 1.8T C ,2.5T C , respectively (when n/ = 3: 
T m in « 1.4T C , 2.0T C , respectively). We note that the curve with /3 = P[2/3] and fi q = in Fig. [Sfb) corresponds to 
the central curve with rif — 3 in Fig. 18 of our previous work [ll[ and to the upper solid- line curve of Fig. 19 of that 
work. 7 

Figs. [7Ja),(b) contain similar results for Ap/T 4 (here only for the value fJb q /T = 0.8). The conclusions about the 
RScl dependence of the results for Ap/T 4 are virtually the same as for p/T 4 . 

It is exactly this independence of v at T down to about 2T C which makes our comparison with lattice data (cf. Fig. [3]) 
much more trustworthy than the simple TPS evaluation. 

In the remaining part we present results for derived quantities, specifically quark number densities and susceptibil- 
ities. 

Fig. [5] contains results for the quark number density (I3"2")) for various values of fi q /T (fii — 0). Included are the 
corresponding results of the lattice calculation of Ref. 21], in the form of points (some with error bars). The values 



7 We mention that a numerical mistake was committed in the mentioned curve of Ref. [Till , in that the power of (j|/ihb) in the program 
there was taken to be five instead of three [cf. Eq. (JTsJ] . The curve with j3 = P[2/3] and [i q = in the present Fig. (6^b) now represents 
the corrected version of the mentioned curve. 
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FIG. 7: As Figs. [6] but for Ap/T 4 instead of p/T 4 . The chemical potential has the values /x^ = 0.8T. 
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FIG. 8: Quark number density n q (divided by T 3 ) as a function of T, at fij = 0, for various values of the ratio fiq/T. Our results are in 
the form of curves. Included as points are the corresponding results of the lattice calculation of Ref. [2]| . 



were obtained by numerical differentiation of our results for p with respect to fi q (with /i/ = and T constant). Again, 
we see that our results in general agree with the lattice results to within 20%, even at low temperatures T ~ T c . 

Finally, in Figs. 0!a),(b) we present the values for the susceptibilities Xq an d Xi [ c f- Eq. ([55)1 ] as a function of 
temperature, for various values of the ratios n q /T, while keeping /ir = 0. The results are presented as various curves, 
and were obtained, at a given T, by numerical evaluation of the double derivatives of our results for the pressure p 
with respect to fi q (at constant /ir = 0) and with respect to \ii (around Mr = 0, at constant fi q ) (24|. In Figs.[HKa),(b) 
we included the corresponding results of the lattice calculation of Ref. [21(, as points with error bars. Our curves, in 
general, give results which are by roughly 20% higher than the lattice results. 

In addition to the aforementioned susceptibilities, the mixed susceptibility Xud which is related to the previous two 
by 

^^ = -4^ (34) 

has been a subject of interest in the literature. Our numerical results give for the above quantity J34|) , at pj = 
(and rif = 2), values of about 4 x 10~ 3 at fj, q = and at temperatures 1 < T/T c < 2 [these values turn negative 
(~ — 10 -3 ) when p q /T = 0.2]. The authors of Ref. [2(| obtained, by their lattice calculations, for the above quantity 
(|34[) at fii = fiq = (and nj = 2) decreasing values as the temperature increases from T c to 1.5T C , and at T/T c = 1.36 
they found a value of 6.6 x 10~ 3 . Our results for this quantity are roughly in agreement with the lattice results of 



Refs. |2Cj, [26J, and with the hard thermal loops (HTL) perturbative estimates of Ref. [35(. The latter estimates give 
for the quantity (|34p . at 1 < T/T c < 2, values between 5 x 10 -3 and 1.4 x 10~ 2 , when using there our values of 
a(2vrT) (with the beta function being BP[2/2]). The lattice quenched results of Refs. [28| give for this quantity values 
~ 10~ 6 , i.e., three orders of magnitude lower. 



15 




U q /T = 1.0 i 

H/T .0.8 

-I : Q 

H/T.0.4 

iyT-0.2 h— *— -I 

ft/T-0.0 i 1 1 



1 1.2 1.4 1.6 1.8 2 1 1.2 1.4 1.6 1.8 2 

T/T c T/T c 
FIG. 9: The susceptibilities (a) \q an d M X/ as a function of temperature, for various values of the ratios [iq/T (/ij = 0). Our results 
are in the form of curves. Included are the corresponding results of the lattice calculation of Ref. [2ll . in the form of points (some with 
error bars). 



Further, we performed numerical calculations of p, Ap and n q in the case of three active flavors nj — 3, with p u — 
= M<? (^-independent) and fx s = 0, with our approach described above. Comparisons with the corresponding lattice 
calculations of Ref. [22j (their Figs. 3 and 6) revealed that, at their values of hb < 0.53 GeV (//b = 3/i g ,riB = n q /3), 
our results for Ap and tib are somewhat higher than theirs, by less than 20% at T/T c > 1.5, and by 30% — 40% at 
T/T c = 1.1. Direct comparison with the lattice results of Ref. ;23f is not possible, as no continuum limit correction 



factor (c^) was applied there. In Ref. [22j, an estimated correction factor c M = 0.446 was applied to Ap. 



IV. SUMMARY 



Within the present paper we extended our recent approach [11] for improving perturbative expressions for the 
quark-gluon pressure (obtained by FTPT) to the case of finite (but still small) quark densities. Thereby, the main 
aim was to find a consistent method for extrapolating the FTPT-based results down to temperatures as low as T c 
(ss 200 MeV). For such low energies, the original TPS's are plagued by huge uncertainties, stemming mainly from 
their strong RScl dependence which itself is partially connected with the occurrence of, at least, two different energy 
scales contributing to the thermodynamic potential under investigation. Therefore, simple FTP series do not permit 
a reliable comparison with existing (low energy) lattice data. Such a check can only be performed if the wild RScl 
dependence is sufficiently tamed. Our method allows such a taming. It rests mainly on two crucial points: Firstly, 
we performed a careful separation of the low-energy from the high-energy contributions to the pressure, which are 
responsible for the (in principle measurable) long- and the short-range behavior, respectively. In this way, we can 
clarify which values of the RScl are the natural ones - they are different for the two parts. Secondly, for each of these 
contributions we identified Pade-related approximants which - besides showing other physically desirable features - 
led to (almost) RScl-stable expressions and thus to predictions which can be safely used down to low temperatures. 
However, the use of the approximants for the low-energy (long-range) contributions at very low temperatures ~ T c is 
only possible if the unphysical perturbative Landau singularities of the QCD coupling parameter at low energies are 
eliminated; we did this by using similar Pade-related approximants for the renormalization-group beta function. As 
a result, we demonstrated that the obtained expressions for the pressure p and the difference Ap — p(vf) — p(0) are 
fairly insensitive to the (as yet) unknown part of contributions of 0(g 6 ) and to variations of the RScl's, both of these 
features being in stark contrast with the TPS expressions. 

Our expressions show a surprisingly good agreement with lattice data - not only for the pressure and its p,f 
dependence but also for derived quantities, in particular susceptibilities. In this context, we note that the unknown 
relative deviations of the low temperature lattice results for p and Ap from the true values are expected to be roughly 
in the range of 10-20%. This is due to the well-known lattice artefacts, in particular the ones connected with finite-/^/ 
effects (truncated Taylor expansion in /z//T). Our Pade-related evaluations give for Ap/T 4 results which are by not 
more than 20% higher than the lattice results when [i q jT < 0.8 and 1 < T/T c < 2, see Fig. [3] 

Our approach is valid only for values of the chemical potentials smaller than the temperature, because only in 
this case dimensional reduction can be applied. Fortunately, present day heavy ion collisions are probing the region 
with values /z/ < 50 MeV, which are small compared to the temperatures T < 5T C typically involved. For other 
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kinematic situations, in particular for small T and larger chemical potentials, different reorganizations of perturbation 
expansions are necessary and have been applied in the literature, the most prominent one being the hard dense loop 
approximation which is genuinely four-dimensional but based on a non-local effective action [29j. Recently, a purely 
diagrammatic calculation of the perturbative QCD pressure (i.e., without involving any effective theory) has been 
performed [3Cj which, at least in principle, should be valid for all kinematic regions. As it should be expected, these 
results - when applied to high temperatures and (relatively) low chemical potentials - are in accordance with those 
of the dimensional reduction approach. 
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APPENDIX A: RELEVANT COEFFICIENTS FOR THE PRESSURE AT FINITE CHEMICAL 

POTENTIAL 

Here we compile expressions for parameters Aj, the and c?e which we obtained from expressions of Ref. @ by 
application of the method of separation of the long-distance from short-distance contributions (i.e., introduction of 
factorization scale Ae: Toe < Ae < 2ttT), as explained in the beginning of Sec. [ill We denote by v the renormalization 
scale (RScl), and g = g{v) in MS scheme. Other notations used in this appendix are: 



2nT 



= — K( fe >(z) = -5>^(z), = 



Uf 



n f 



In' £(— n) 



C(-n) 
C(-n) 



(Al) 
(A2) 

(A3) 



where H(z) and H(£, z) are the aleph functions defined in Refs. Q via digamma functions and derivatives of the 
Riemann zeta functions. 

Coefficients Aj (j = 1, ...7) and Aj , which appear in equations Eqs. I®-© and (f5])- (fTU|) . and later in Eqs. ([15]) - 
(I20p . are obtained from the "matching parameters" ag; (i — 1, ...7) in Ref. [6j by separating appropriately the parts 
proportional to hiv c (In A in Ref. [6j), or proportional to 1/e, from the remaining (y c - and e-independent) parts. They 
take the form 



A x -- 
A 2 = 
A 4 -- 
A 5 = 

+ 
A 7 -- 



7T 

45 
1 

_ 6 



3nf ( - + 30/i 2 + 60/i4 



1 + — n f (5 + 72^2 + 144^) 



1 + ~n f (1 + 12/12) 



2 (-In 2 + ln'C(-l)) + ^n f (l - 21n2)(l + 12/I 2 ) + 4n / H^(l, z) , 
j ^ [90 - 396 In 2 + 66 7E (6 + n f ) + n f (3 - 42 In 2) + 2nj (1 + 2 In 2)] 

n/^2 I 6 + 44 7E - 44 In 2 + n f -(l + 2 In 2) ] + -n f (n f + 6)H (0) (z) + n 2 f -H (2) (z) 
V 3/9 1 3 



22( 7E - In 2) +l+n / -rn2 + n / ^ ( °)(^) 



(A4) 
(A5) 
(A6) 
(A7) 

(A8) 
(A9) 
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.4 



2A 4 , 



4" } = 22 + -n f - -n 2 + 4(11 - -n f )n f Jl 2 , 



.4 



22 



(A10) 
(All) 

(A12) 



The most complicated coefficient is A3, which emerges both in Eq. ([8]) and later in Eq. (|20|) via the quantity k 
[Eq. (|2"T])]. It can be expressed as 



where 



A 3 ,i 



18 



H6 „ 38, tJ ., . 220, , A . „, 
— + 4 7 E - -7T- In C -3 + — In' ( -1 
5 3 3 



^ + 8 7E + 2(127 + 48 7E )M2 - 644/i 4 + ^ In' £(-3) + -(11 + 156Jh) In' CM) 
60 15 3 



+24 (52K (0) (3, z) + 144iN (1) (2, z) + 17K (0) (1, z) - 92K (2) (1, z) + 4iK (1) (0, z) 



-(35 + 4727I 2 + 1328j5 4 ) - 24(1 - 4/l 2 ) In' f (-1) 



-144 (l2iK (1) (2, 2) - 2N (0) (1, z) - 16tt (2) (l, z) - iN (1) (0, z) - i4N (3) (0, z)) 



1 64 32 

- + 4 7E + 8(7 + 12 7E )M2 + 112ft - — In' C(-3) - — (1 + 12/i 2 ) In' C(-l) 

3 15 3 



-96 



(8H( 0) (3, z) + il2N«(2, z) - 2K(°)(1, z) - 4:^(1, z) - iN«(0, z)) 



(A13) 



(A14) 



and 



-4, 



3.2 



4^ 2(1 + 7E )7i^ - K(3, */ + z g ) - H(3, z/ + z g *) - 1471/ (K(2, z f + z g ) + K(2, z f + z* g )) 
f,g L 

+477 2 N(1, z f ) + (77/ + 77 g ) 2 K(l, z f + z g ) + (77/ - 77 g ) 2 N(l, z f + z* g ) + i477 2 77/K(0, z f ) , 



A 



97 . 169 5 2 \ . „ /„„ . 4 \ . „ / . . . 8 



54 



+ "36~ n/ + ) + mUf ' 50 + ) + M4n/ ' _44 + * Uf 



The constants ac, ctMi and «M2 were obtained in Ref. 



43 157 



43 491 



96 6144 



^ « 0.195715, OMi = — - — — 7T Z « 0.555017, a M2 = -- , 

32 6144 3 



and /3m 1 in Ref. [311 ] : 



Ml 



-1. 391512 



(A15) 



(A16) 



(A17) 



(A18) 



APPENDIX B: BETA FUNCTIONS OF THE BOREL-PADE TYPE 



In this Appendix we present various resummations of the QCD f3 functions as functions of x = a(Q 2 ) = a s (Q 2 )/7r, 
in the MS scheme. Further, the corresponding running of a(y 2 ) as function of x = v 2 (in GeV 2 ) is given, in the 
various cases, always normalized to a(m 2 ) = 0.334/7T = 0.106316, cf. Eq. (|29|) . We will denote the squared RScl v 2 
here as Q 2 (= —q 2 > 0) to emphasize the space-like character of the corresponding four- vector q. 
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ft = 57 H; t*-»/ + xt»/ > ( B3 ) 



The four-loop renormalization group equation (RGE) for a(Q 2 ) = a s (Q 2 )/ir is: 

g2 d^) = _ ftfl2(Q2) [x + cifl(Q2) + C2fl2(g2) + C3a 3 (Q2)] ) (B1) 

where c,- = Pj/0o (j > 1)- The one- and two-loop coefficients /3o and /3j [12, 1331 are scheme independent; in MS 
scheme [34j the three- and four- loop coefficients /3 2 and /3 3 were obtained in Refs. |35l. l36j , respectively 

A> = ^ (n - fn/) , A = ^ (l02 - yn/) , (B2) 

1 /2857 5033 325 2 
64 {~2 W Uf + IT"' 

fa = -^-(29243.0 - 6946.30 n f + 405.089 n 2 f + 1.49931 nl) , (B4) 
256 

and n/ is the active number of quark flavors. 

The solution of the RGE (|B1[) . at low Euclidean energies Q (with either ny=3or?iy = 2) has the known 
unphysical Landau singularities, i.e., singularities of a(Q 2 ) for Q 2 < Qp oic - For example, if we choose the realistic 
value of a s {Q 2 = to 2 ) = 0.334, the singularity in a(Q 2 ) appears already at Q 2 ok , ~ 0.66 2 GeV 2 « 0.44 GeV 2 when 
n f = 3, and Q 2 olc w 0.75 2 GeV 2 w 0.57 GeV 2 when n/ = 2 - cf. Figs. EUa) andE^a). The main reason for this 
unphysical behavior of a(Q 2 ) is the truncated perturbation series (TPS) form of the beta function j3(x = a) - the 
right-hand side of RGE (|B1[) . Such a form of (3(x) has its origin in the perturbative approach (powers of x = a). The 
TPS |/3(x)| grows out of control when x — a increases - cf. Figs. ITUT a) and [T4T a). This leads to the appearance of the 
nonphysical singularities in a(Q 2 ) at low positive Q 2 < Qp olc - 

These singularities prevent us from using, in the traditional perturbative QCD (pQCD), the coupling at squared 
energies Q 2 < Q 2 i c . The problem can be avoided by certain resummations of the TPS (3 function, i.e, by finding 
such a /3(a) function whose Taylor expansion around a = up to ~ a 5 reproduces the TPS /3(a) of Eq. (jBlj) and, at 
the same time, /3(a) remains more under control when a increases. 

One possibility is to construct diagonal or near-to-diagonal Pade approximants based on the TPS /3(a). For example, 
the Pade 

a 2 

P[2/3]fl(o) = -Pot- r-= ; — ~ — r-^r (B5) 

L ' JPW [l-cia+(c 2 -c 2 )a 2 + (-c? + 2cic 2 -c 3 )a 3 ] v ' 

gives us an expression which, up to o k 0.3, behaves well - cf. Figs. [TUT b) and [FiTb) . However, around a « 0.3, this 
/3(a) goes abruptly out of control, because the Pade expression has a pole there. The corresponding running coupling 
a(x = Q 2 ) achieves singularity already at Q 2 olc w 0.81 2 GeV 2 w 0.65 GeV 2 for n f = 3, and Q 2 olo w 0.92 2 GeV 2 w 
0.85 GeV 2 for nj = 2 - cf. Figs. IT2Tb) and fTBT b). In contrast to the TPS /3(a) case, however, a(Q 2 ) seems to be well 
under control now for virtually all Q 2 larger than Q 2 olc . 

Another possibility, which avoids the aforementioned pole problem of the Pade fi(x = a), would go in the direction 
of resumming first the Borel transform ~Bp(y) (the latter has in general significantly weaker singularities than /3), 
and then applying the inverse transformation via a Borel integration. For example, we can try to apply diagonal or 
close-to-diagonal Pade resummation to B^(y). The Borel transform is 

Mv) = ~ /3o (n + Cl |7 + C2 |T + C3 |y + ---) > ( B6 ) 

and the Pade P[2/2] and P[l/3] resummations of the above TPS are 

P[2/2] B (y) = -/3o 1 J / . + r f. 2 2 , (B7) 
1 + hy + t 2 y 2 

P[l/3] B (iO = -/3o— J <a , 3 , (B8) 

1 + siy + s 2 y z + s 3 y 6 

where the coefficients rj,tj, Sj are unique functions of Cfe's such that reexpansion of (|B7[) and (|B8[) reproduces expansion 
(|B6p up to (and including) ~ y 4 term: 

r 2 = (l/2)(3c?-4c lC2 + c 3 )/e , 

h = (l/2)(-2c lC2 + c 3 )/e , ta - (l/12)(4c 2 -3c lC3 )/cf , (B9) 
Sl = (l/2)ci , s 2 - (l/12)(3c 2 - 2c 2 ) , s 3 = (l/24)(-3c 3 + 4 Cl c 2 - c 3 ) , (B10) 
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and we used the notation £ = (3c 2 — Ic-i). However, the inverse Borel transformation 



(3{x) = dy exp(-y/x)Bp(y) 
Jo 



(Bll) 



cannot be constructed by inserting here directly the Pade expressions (|B7|) or (|B8|) for the integrand. This is so 
because the latter expressions have poles on the positive axis: P[2/2]b(?/) at y p . ps 1.13, 1.01 for = 3, 2, respectively; 
P[1/3]b(2/) a t Up. ~ 0.94,0.88 for rif — 3,2, respectively. These "infrared renormalon" singularities imply ambiguities 
in the integration: 5f3(a) ~ exp(— y p /a) [~ S(Q 2 da(Q 2 ) / dQ 2 ) ~ <5a(Q 2 )]. This implies an ambiguity in the coupling 
a{Q 2 ): 5a(Q 2 ) ~ (A^/Q 2 )^- = (A^/Q 2 )''. Numerically, for P[2/2] B , V ~ 2.55,2.44 for n f = 3,2, respectively; 
for P[1/3]b, f) ~ 2.11,2.13 for n/ = 3,2, respectively. We can fix the above ambiguity by choosing a specific recipe 
for the Borel integration over the pole y p .. We will choose the Principal Value (PV) prescription 



poo±ie 

BP [*/j]/3(20 = Re / dy exp( 

J±ie 



-y/x)P[i/j] B (y) 



(B12) 



where [i/j] = [2/2] or [1/3]. Numerically, this is difficult to implement, as e — > +0 and we approach the pole y p . 
down to the distance s during the integration. However, we can use the Cauchy theorem, and the fact that the Borel 
transforms (|T37|) and (|B8|) do not have any poles in the complex semiplane Re(y) > except the aforementioned 
y p . > 0. This allows us to avoid the vicinity of the pole, for example by integrating along a ray y = r exp(— where 
(f> is any small but finite positive fixed angle (cf. Refs. [xTL fl9|) 



BP[i/i]/3(a:) = Re|exp(-i0) / dr exp(-y/x)P[i/j] B (y) 



y=r ex.p(—i4 



(B13) 



This approach, which is numerically stable, gives us for the Borel-Pade(BP)-resummed /3(x)-functions values which 
are surprisingly non-singular and achieve at x = a rs 1. value zero (infrared fixed point) - cf. Figs. [11] and 1151 
Integration of the RGE with these BP-rcsummed four-loop MS j3 functions, with the phenomenologically acceptable 
initial condition a s (m 2 ) = 0.334 [Eq. (|2"9"|) ]. gives us a(x = Q 2 ) running couplings which are presented in Figs. [TBI and 
[T71 Both choices BP [2/2] and BP [1/3] give similar behavior for a(Q 2 ). There are no Landau singularities present 
any more, and the coupling is analytic in the sense that there are no unphysical singularities on the space-like axis of 
the squared momenta q 2 = —Q 2 . Furthermore, due to the aforementioned zero of the BP (3 functions, the coupling 
a{Q 2 ) remains finite down to Q 2 — where it has a value m 1. The obtained "analytized" coupling a(Q 2 ) probably 
represents a version of analytic QCD. Therefore, the skeleton-motivated method of Refs. [37j can be applied as a 
alternative way of evaluating the low-energy QCD observables. 
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FIG. 10: MS beta functions -/3(x)//3 , for n f = 3, whose TPS is TPS(x) = 1 + cix + c 2 x 2 + c 3 x 3 : (a) TPS(x), and (b) Pade P[2/3](x) 
here stands for a(Q 2 ). 
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FIG. 12: MS running of a(x = Q 2 ), for n f = 3, when the beta functions arc (a) TPS, and (b) Pade P[2/3]. 



(beta=BP [2/2] ,bMS, nf=3) 




a f x ] (b) (beta=BP [1/3] ,bMS,nf=3) 




0.2 0.4 0.6 0. 



0.2 0.4 0.6 0. 



-0.2 L -0.2 L 

FIG. 13: MS running of a(x = Q 2 ), for n f = 3, when the beta functions arc (a) Borcl-Pade BP[2/2], and (b) BP[l/3]. 
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FIG. 14: MS beta functions -/3(x)//3 , for n f = 2, whose TPS is TPS(x) = 1 + cix + c 2 x 2 + c 3 x 3 : (a) TPS(x), and (b) Pade P[2/3](x) 
here stands for a(Q 2 ). 
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